V us from kaon decays 
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Recent measurements of kaon decays contributing to the determination of \V US \ are summarized, 
and up-to-date evaluations of |V us |/+(0) and \V US \ are presented. 



I. INTRODUCTION 



II. K L DECAY RATE DATA 
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At present, the first-row constraint |V uc i| 2 + | V^, s | 2 + 
\V u b\ 2 (with \ V u b\ 2 negligible) offers the most precise test 
of CKM unitarity. Up until 2002 (and for the 2004 PDG 
evaluation [l|), the evaluation of \V US \ from older Ktz 
data gave a 2.3c hint of unitarity violation in the first- 
row test. The 2003 measurement of BR(K+ 3 ) by BNL 
E865 gave a value for \V US \ consistent with unitar- 
ity. In the period 2004-2006, many new measurements 
of BRs, lifetimes, and form-factor slopes were announced 
by KLOE, KTeV, ISTRA+, and NA48. All of these new 
measurements are distinguished from the older measure- 
ments in that they are based on much higher statistics, 
and in that radiative corrections are applied consistently. 
Late-2005 evaluations of \ V US | by the CKM working group 
on the first row [|| and in the 2006 PDG review [i[ include 
many, but not all, of these important developments. We 
present an up-to-date evaluation that includes prelimi- 
nary results presented at this conference. 

\V US \ is evaluated from Kg 3 data using the relation 
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where Ck is a Clebsch-Gordan coefficient; Mk is the 
appropriate kaon mass; Sew = 1-0232 is the universal 
short-distance electroweak correction; /+(0) = f+ * (0) 
is the value of the hadronic matrix element (form factor) 
for the K° — > ir~ transition at zero momentum transfer 
to the leptonic system; Iki is the phase-space integral of 
the normalized form factor, which depends on the values 



of one or more slope parameters A; A^ 7 ^ and A^-" are 
5?7(2)-breaking and long-distance electromagnetic cor- 
rections; and the subscripts K and £ indicate dependence 
on the kaon charge and lepton flavor. Our main result 
is a current average value for |U us |/+(0). We first review 
the recent data on Kn 3 rates (BRs and lifetimes) and 
form-factor slopes. 
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Numerous measurements of the principal Kl BRs, or 
of various ratios of these BRs, have been published re- 
cently. For the purposes of evaluating (0), these 
data can be used in a PDG-like fit to the Kl BRs and 
lifetime, so all such measurements are interesting. 

KTeV has measured five ratios of the six main Kr, 
BRs @ . The six channels involved account for more than 
99.9% of the Kl width and KTeV combines the five mea- 
sured ratios to extract the six BRs. We use the five mea- 
sured ratios in our analysis: BR(K^/ K e3 ) = 0.6640(26), 
BR(7r+7r-7r°/A' e 3) = 0.3078(18), BR{ir+n- / K e3 ) = 
0.004856(28), BR(3ir°/K e3 ) = 0.4782(55), and 
BR(27r°/37r°) = 0.004446(25). The errors on these 
measurements are correlated; this is taken into account 
in our fit. 

NA48 has measured the ratio of the BR for K e3 decays 
to the sum of BRs for all decays to two tracks, giving 
BR(A- e3 )/(l - BR(3tt )) = 0.4978(35) @. From a sep- 
arate measurement of BR(K L 3n°)/BR(K s -> 2tt°), 
NA48 obtains BR(37r°)/r L = 3.795(58) MHz [|. 

Using cj) — > KlKs decays in which the K$ decays to 
7r + 7r~, providing normalization, KLOE has directly mea- 
sured the BRs for the four main Kl decay channels 
The errors on the KLOE BR values are dominated by 
the uncertainty on the Kl lifetime tl] since the depen- 
dence of the geometrical efficiency on tl is known, KLOE 
can solve for tl by imposing J2 X BR(A^l — ► x) = 1 (using 
previous averages for the minor BRs) , thereby greatly re- 
ducing the uncertainties on the BR values obtained. Our 
fit makes use of the KLOE BR values before applica- 
tion of this constraint: BR(K e3 ) = 0.4049(21), BR(A^ 3 ) 
= 0.2726(16), BR(K e3 ) = 0.2018(24), and BR(K e3 ) = 
0.1276(15). The dependence of these values on tl and 
the correlations between the errors 0] are taken into ac- 
count. 

KLOE has also measured tl directly, by fitting the 
proper decay time distribution for Kl — > 37r° events, for 
which the reconstruction efficiency is high and uniform 
over a fiducial volume of ^0.4Al. They obtain tl = 
50.92(30) ns [l|. 

There are also two recent measurements of 
BR(7r + 7r~/A^3), in addition to the KTeV measure- 
ment of BR(tt + it~ /K e3 ) discussed above. KLOE 
obtains BR(7r+7r-/if M3 ) = 7.275(68) x 10~ 3 11], while 
NA48 obtains BR(tt + it- /K e3 ) = 4.826(27) x 10~ 3 
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TABLE I: Results of fit to K L BRs and lifetime 



Parameter Value S 

BR(Kes) 0.40563(74) LT 

BR(_ftT M 3) 0.27047(71) 1.1 

BR(3tt°) 0.19507(86) 1.2 

BR(7T + 7r"7r°) 0.12542(57) 1.1 

BR(7r + 7r-) 1.9966(67) x 10" 3 1.1 

BR(2tt°) 8.644(42) x 10" 4 1.3 

BR( 77 ) 5.470(40) x 10" 4 1.1 

tl 51.173(200) ns 1.1 



[12j |. All measurements are fully inclusive of inner 
bremsstrahlung. The KLOE measurement is fully 
inclusive of the direct-emission (DE) component, DE 
contributes negligibly to the KTeV measurement, and a 
residual DE contribution of 0.19% has been subtracted 
from the NA48 value to obtain the number quoted 
above. For consistency, in our fit, a DE contribution of 
1.52(7)% is added to the KTeV and NA48 values. Our 
fit result for BR(7r + 7r~) is then understood to be DE 
inclusive. 

In addition to the 14 recent measurements listed above, 
our fit for the seven largest Kl BRs and lifetime uses four 
of the remaining five inputs to the 2006 PDG fit and the 
constraint that the seven BRs sum to unity. The results 
are summarized in Table [TJ 

Our fit gives x 2 /ndf = 20.2/11 (4.3%), while the 2006 
PDG fit gives x 2 /ndf = 14.8/10 (14.0%). The differences 
between the output values from our fit and the 2006 PDG 
fit are minor. The poorer value of x 2 /ndf for our fit 
can be traced to contrast between the KLOE value for 
BR(37T°) and the other inputs involving BR(3"7r°) and 
BR(7T°7r )— in particular, the PDG ETAFIT value for 
BR(7r 7r°/7r+7r"). The treatment of the correlated KTeV 
and KLOE measurements in the 2006 PDG fit gives rise 
to large scale factors for BK(K e3 ) and BR(37r°); in our 
fit, the scale factors are more uniform. As a result, our 
value for BB.(K e3 ) has a significantly smaller uncertainty 
than does the 2006 PDG value. 



III. K s DECAY RATE DATA 

Two KLOE measurements provide the BR for the K e3 
decay of the K$ with enough precision to be of inter- 
est for the determination of \V US \. Using <f) — * K^Kg 
decays in which the Kl is recognized by its interac- 
tion in the experiment's EM calorimeter, KLOE mea- 
sures BR(X e3 /7r+7r-) = 10.19(13) x 10~ 4 [3. KLOE 
has also recently obtained BR(7r+7r"/7r 7r ) = 2.2459(54) 
[Til ], where this value is an average including the ex- 
periment's 2002 result. These measurements completely 
determine the main K$ BRs, and in particular give 
BR(K e3 ) = 7.046(91) x 10" 4 [l|. Our evaluation of 
\V us \f+(0) uses this value, together with the lifetime 
value t s = 0.08958(5) ns from the PDG fit to CP pa- 



TABLE II: Results of fit to K± BRs and lifetime 



Parameter Value S 

BR(A- m2 ) 63.442(145)% O" 

BR(7T7r°) 20.701(108)% 1.3 

BR(7T7r7r) 5.5921(305)% 1.0 

BR(A" e3 ) 5.121(38)% 1.6 

BR(ie M 3) 3.3855(203)% 1.2 

BR(7T7r 7r°) 1.7592(234)% 1.1 

r± 12.3840(213) ns 1.8 



rameters, which is highly constrained by measurements 
from NA48 [H and KTeV [H. 



IV. A'± DECAY RATE DATA 

There are several new results providing information on 
rates. These results are mainly preliminary and have 
not been included in previous averages. 

NA48/2 has recently submitted for publication 
measurements of the three ratios BR(if e3 /7T7r°), 
BR^s/tttt ), and BR(K ^ / K e3 ) [13, d- These mea- 
surements are not independent; in our fit, we use the 
values BR^s/tttt ) = 0.2496(10) and BR{K^/ttt: q ) = 
0.1637(7) and take their correlation into account. 

ISTRA+ has also updated its preliminary value for 
BR(i^ e3 /7r7r°). They now quote BR(^ e3 /7r7r°) = 
0.2449(16), as reported at this conference jlQj ] - 

KLOE has measured the absolute BRs for the K e3 and 
decays [2(| ■ In 4> ~ * K + K~ events, K + decays 
into ijlv or irir are used to tag a K ~ beam, and vice 
versa. KLOE performs four separate measurements for 
each K13 BR, corresponding to the different combinations 
of kaon charge and tagging decay. The final averages 
are BR(K e3 ) = 5.047(92)% and BR(X m3 ) = 3.310(81)%. 
Our fit takes into account the correlation between these 
values, as well as their dependence on the K ± lifetime 

The world average value for t± is nominally quite pre- 
cise; the 2006 PDG quotes r± = 12.385(25) ns. How- 
ever, the error is scaled by 2.1; the confidence level for 
the average is 0.2%. It is important to confirm the 
value of t±. A preliminary measurement from KLOE, 
t± = 12.367(78) ns [9j, agrees with the PDG average, al- 
though at present the KLOE uncertainty is significantly 
larger. 

Our fit for the six largest K BRs and lifetime makes 
use of the results cited above, plus the data used in the 
2006 PDG fit, for a total of 30 measurements. The six 
BRs are constrained to sum to unity. The results are 
summarized in Table HT1 

The fit quality is poor, with x 2 /ndf = 49/24 (0.21%). 
However, when the five older measurements of t± are re- 
placed by their PDG average with scaled error, x 2 / n df 
improves to 31.3/20 (5.1%), with no significant changes 
in the results. Tension between the new measurements 
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FIG. 1: Evolution of average values for main K± BRs. 



involving BR(A" e 3) and the older measurements is partly 
responsible for the poor fit quality (note the scale factor 
of 1.6 on BR(if e 3) from our fit). Because of uncertain- 
ties concerning the treatment of radiative corrections, it 
would be desirable to eliminate some of the older mea- 
surements, as both we and PDG do for the Kl fit. How- 
ever, for K , all of the new measurements involve the 
BRs or the ratios BR(X£3/7T7r°). This leads to large cor- 
relations between the BRs and BR(7r7r°), and makes 
the fit unstable when the older data are excluded. Both 
the significant evolution of the average values of the K(z 
BRs and the effect of the correlations with BR(7T7r°) are 
evident in Fig. [TJ 



V. FORM-FACTOR SLOPES 

For K e s decays, recent measurements of the quadratic 
slope parameters of the vector form factor (A' + , A" ) are 
available from KTeV [H H2, KLOE 23], ISTRA+ 
[lH, and NA48 [25[. For decays, results of fits 
using the quadratic parameterization for the vector form 
factor (A' + ,A+) and the linear parameterization for the 
scalar form factor (Ao) are available from KTeV [2l[, IS- 
TRA+ [H, HI, and NA48 [jj, Hz|. ISTRA+ measures 
K~ decays; the other experiments measure Kl decays. 

The K e3 data are summarized in Fig. [2] The uncer- 
tainties on the values of A' + and A" reported by each 
experiment are highly correlated. This is an intrinsic 
property of the quadratic parameterization. The form- 
factor slopes represent a small modification of the kine- 
matic density of the matrix element for K — > 7r transi- 
tions; in addition, sensitivity to A" is poor because the 
kinematic density of the matrix element drops to zero 
at large t, where the form factor itself is maximal. Tak- 
ing the correlations into account, we obtain a good fit 
to the data in Fig. 03 A' + = 24.15(0.87) x 1CT 3 and 
A" = 1.57(0.38) x 10~ 3 , with p{\'+,\'+) = -0.941 and 
X /ndf = 5.3/6 (51%). The significance of the quadratic 
term is greater than 4cr. The fit result is shown as the 
yellow ellipse in Fig. [2] 

The K^3 data are summarized in Fig. [31 which also 
shows the results for A+ and A" from our K e3 fit, and the 
results of our combined fits to K e3 and K^ 3 data. While 
the results for A^ and A" from K^ 3 decays are fairly con- 
sistent and reasonably agree with those from K e3 decays 
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FIG. 2: Recent measurements of form-factor slopes for K e 3 
decays, with fit. 



(Fig. [21 left panel), the results for Ao show poor consis- 
tency (right panels). The most discrepant measurement 
is that of NA48. In both the A+-Ao and A" -Ao views, a 
locus of values within the la confidence contours from 
ISTRA+ and our K e3 fit lies along the major axis of 
the KTeV confidence ellipse. As a result, a fit to all data 
(K e3 and K^ 3 ) excluding only the results from NA48 
gives x 2 / ndf = H-9/9 (21.7%). The results of this fit are 
indicated by the cyan ellipses in Fig. [3] When the NA48 
results are included, x 2 / n df increases to 52/12, for 
a probability of less than 10 -6 . The results are shown as 
the yellow ellipses in Fig. [3) 

Since at the moment there is no a priori reason to 
exclude the NA48 results, we base our estimate of 
|Ki S |/+(0) on the fit to all data. However, we scale the 
errors on the slope parameters to reflect the inconsistency 
in the input data set. We obtain A^_ = 24.84(1.10) x 10~ 3 
(S = 1.4), A" = 1.61(0.45) x 10" 3 (S = 1.3), and 
Ao = 13.30(1.35) x 10~ 3 (S = 2.1), with p(\'+,\+) = 
-0.944, p(A' + ,A ) = +0.314, and p(A£,A ) = -0.420. 
From these results, we calculate the phase-space inte- 
grals for use with Eq. CQ) to be I(K%) = 0.15452(29), 



) = 0.15887(30), 
0.10501(35). 



H3J 



0.10207(34), and 
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FIG. 3: Recent measurements of form-factor slopes for Kts decays, with fits. 
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TABLE III: Summary of 5"(7(2)-breaking and EM corrections 



Mode 


A SU(2) 


a KM 


K% 




+2.31(22)% 


+2.31(22)% 


+0.52(10)% 
+0.03(10)% 
+0.95(15)% 
+0.2(4)% 



VI. EVALUATION OF \V us \f+(0) 



TABLE IV: Values of |V us |/+(0) from data for different decay 
modes, with error breakdown 









App 


rox contrib to % 


err 


Mode 


\Vus\f+(0) 


% err 


BR 


T 


A 


i 


K L e3 


0.21639(55) 


0.25 


0.09 


0.19 


0.10 


0.09 


Ktfj3 


0.21649(68) 


0.31 


0.10 


0.18 


0.15 


0.17 


K s e3 


0.21555(142) 


0.66 


0.65 


0.03 


0.10 


0.09 


K ± e3 


0.21844(101) 


0.46 


0.38 


0.11 


0.24 


0.09 


K ± /i3 


0.21809(125) 


0.57 


0.31 


0.10 


0.45 


0.17 



The strong 5*[/(2)-breaking and EM corrections to 
/+(0) used in our evaluation of |V^ s |/ + (0) are summa- 
rized in Table QIIJ For A su ^ , we use the standard 
value from Ref. |28|. All of the recent Ki§ BR measure- 
ments are fully inclusive of final-state radiation; the val- 
ues we use for A™ were calculated for the fully-inclusive 
rates. Specifically, for K° 3 and Kf 3 , the values were ob- 
tained from the chiral perturbation theory calculation 
of Ref. [H, updated using the evaluations of the low- 
energy constants from Ref. H{| For K® 3 , the value is 
from Ref. and was obtained using a generator imple- 
menting a hadronic- model calculation. For K^ 3 , we do 
not know of any complete estimate. Our value is loosely 
based on Ref. l3ll . with a very generous error estimate. 
However, with this treatment, A™ gives the largest sin- 
gle contribution to the uncertainty on |V^ s |/+(0) from 
this mode. 

We use the results of the fits and averages described 
in Sees. HlHVl together with the corrections in Table HTT1 
to evaluate the quantity |V^ s |/+(0) for each of the five 
Ki3 decay modes as listed in Table IIVI The table also 
gives the fractional uncertainty for each determination 
of |Vu S |/+(0), as well as the approximate breakdown of 



the contributions to the uncertainty from the various in- 
puts (these contributions are understood to be added 
in quadrature). The source of the limiting uncertainty 
varies from mode to mode. 

The average value of |T4s|/+(0) from all modes, as ob- 
tained from a fit with correlations taken into account, is 
0.21673(46). The fit gives % 2 /ndf = 4.2/4 (38%). 

Although the value of y 2 /ii&i from the fit is satis- 
factory, one notes that the values of \V us \f + (0) for the 
two charged modes seem to be higher than those for 
the neutral modes. To quantify this, we perform the 
fit separately for charged and netural modes, using the 
results of the overall fit to the form-factor slope data 
of Sec. El and the S'(7(2)-breaking correction from Ta- 

SU(2) 



bleUm A^^ = 2.31(22)%, for the charged modes. We 
obtain |V A us |/+(0) = 0.21635(50) for neutral modes and 
0.21832(94) for charged modes, a 1.9a difference. If we 
take this as a suggestion that the SU (2)-breaking correc- 
tion may be underestimated, and perform the fit leaving 
free the value of \V US |/+ (0) for neutral modes and an em- 
pirical value of A su( > 2 \ we obtain Afip (2) = 3.24(43)%. 
Comparison of the values of |V^ s |/+(0) for K e ^ and 
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FIG. 4: Results of fits to |Kd|, \V U s\, and |I4 s |/|Kd|. 

modes provides a test of lepton universality. Specif- 
ically, 

_ [|K a |/+(0)]^ 3 ,cxp _ I> Ie3 (1 + S e3 ) 

r " e [|K s |/+(o)]2 3 ,e XP r e3 '/ M3 (i + M' 

By comparison with Eq. {!]), r Me is equal to the ratio 
wrtn 9i the coupling strength at the W ^> Iv ver- 
tex. In the Standard Model, r^ e = 1. We calculate r Me 
separately for charged and neutral modes. The results 
are compatible; the average value is r^ e = 1.0003(53). 
We note that the sensitivity of this test of lepton univer- 
sality is beginning to approach that obtained in 7r — > Iv 
decays: (r eM ) T(!2 = 0.9966(30), as reviewed in Ref . l32l. 

VII. EVALUATION OF |V US | 

To test CKM unitarity, a value for /+(0) is needed. 
In the chiral expansion, /+(0) = 1 + f p * + f p e + • • • . 
A finite one- loop contribution gives rise to / p 4j its cal- 
culation [H, |34| involves essentially no uncertainty. A 
quark- model estimate of the contribution at 0(p 6 ) was 



originally performed by Leutwyler and Roos 33j. More 
recently, a complete 0(p 6 ) calculation in chiral pertur- 
bation theory showed that meson-loop contributions are 
sizeable [3|| . Analytic estimates of /+ (0) were surveyed 
at this conference (3||, and lattice evaluations of f+(0) 
are rapidly improving in precision (see, e.g., Ref. 1371 ). 
Until a definite consensus emerges on a new reference 
value for /+(0), however, we use the original estimate 
of Leutwyler and Roos, /+(0) = 0.961(8), which is sup- 
ported by lattice calculations [13] • This gives \V US \ = 
0.2255(19). 

Using the recent evaluation of \V U d\ from + — > + 
nuclear beta decays (Ref. [H; see also Ref. |H), \V u d\ = 
0.97377(27), one has \V ud \ 2 + \V US \ 2 - 1 = -0.0009(10), 
a result compatible with unitarity. 

Marciano [40J has observed that T^K^/T^^) can be 
precisely related to the product (\V us \/\V ud \) 2 (fK/U) 2 - 
The recent measurement BR(A" + — » n + v) = 0.6366(17) 
from KLOE [4l|, together with the preliminary lattice 
result f K /U = 1.208(2)(±J 4 ) from the MILC Collabora- 
tion 0, gives |V ua |/|Kd| = 0.2286(1^). This ratio can 
be used in a fit together with the values of \V u d\ from 
Ref. [H and \V US \ from decays as above. Using the 
value for \V ua \ obtained with /+(0) = 0.961(8), the fit 
gives \V ud \ = 0.97377(27) and \V US \ = 0.2245(16), with 
% 2 /ndf = 0.77/1 (38%). The unitarity constraint can 
also be included, in which case the fit gives % 2 /ndf = 
3.10/2 (21.2%, or 1.2a). Both results are illustrated in 
Fig. [4] 

In summary, from data we obtain |V^ s |/+(0) = 
0.21673(46), where the uncertainty amounts to 0.21%. 
The dominant contribution to the uncertainty on \ V US \ is 
from /+(0). Whether \V US \ is evaluated from Ki 3 data 
alone or with the additional constraint from decays, 
the first-row unitarity test is satisfied at about the la 
level. 
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